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I . Abstract. The interaction of a weakly nonlinear Alfvenic pulse with an Alfven speed 
" ' inhoniogeneity in the direction perpendicular to the magnetic field is investigated. Iden- 
^~^ tically to the phase mixing experienced by a harmonic Alfven wave, sharp transverse 
^^ gradients are generated in the pulse by the inhomogeneity. In the initial stage of the 

00 : 

iy~j ,. evolution of an initially plane Alfvenic pulse, the transverse gradients efficiently gener- 

,«-^ • ate transversely propagating fast magnetoacoustic waves. However, high resolution full 

"^ ' MHD numerical simulations of the developed stage of the pulse evolution show that the 

^H ' generation saturates due to destructive wave interference. It is shown that the weakly 

"T"^, non-linear description of the generated fast magnetoacoustic wave is well described by 

^ the driven wave equation proposed in Nakariakov, Roberts & Murawski (1997), and a 
-*— > ■ 

^ ■ simple numerical code (2D MacCromack) , which solves it with minimal CPU resources 

CZ ■ 

^ I produces identical results to those obtained from the full MHD code {Lare2d, Arber et 

•I— I . X , . 

k^ , al. 2001). A parametric study of the phenomenon is undertaken, showing that, contrary 
i3 : to one's expectations, steeper inhomogenities of the Alfven speed do not produce higher 
saturation levels of the fast wave generation. There is a certain optimal gradient of the 
inhomogeneity that ensures the maximal efficiency of the fast wave generation. 
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1. Introduction 

Many astrophysical objects and space plasma systems are significantly influenced by 
magnetohydrodynamic (MHD) waves. It is enough to mention that MHD waves consti- 
tute the basis of many current theories of solar and stellar wind acceleration, the heating 
of solar and stellar coronae and the stability of molecular clouds, etc. In solar physics 
interest in MHD waves has increased significantly with the successful detection of MHD 
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waves in coronal loops and plumes by recent observations from the SOHO and TRACE 
spacecrafts. This, in turn, has generated lively interest in MHD waves in the theoret- 
ical solar physics community in the context of coronal plasma heating and solar wind 
acceleration (see Roberts 2000 for a review). 

A significant aspect in MHD wave dynamics, which has attracted considerable at- 
tention recently is the interaction of waves with plasma inhomogeneities. Because of the 
essentially anisotropic nature of MHD waves, they are affected by the inhomogeneities in 
many ways. In particular, structuring of the medium leads to the appearance of long wave 
length dispersion, the existence of local resonances associated with e.g. resonant absorp- 
tion, and Alfven wave phase mixing (see, e.g. Roberts 1991 and Roberts & Ulmschneider 
1997). 

In a majority of astrophysical applications, almost incompressible Alfven waves are 
subject to very weak dissipation. Therefore, physical mechanisms which can create small 
scales in the waves, enhancing the dissipation, are very important. One of these mecha- 
nisms is phase mixing of Alfven waves interacting with a plasma inhomogeneity. When 
an initially plane Alfven wave propagates along a straight magnetic field, the presence of 
an Alfven speed inhomogeneity across the field leads to distortion of the wave front. In 
other words, each magnetic interface supports Alfven waves propagating with its own lo- 
cal Alfven speed. After a while, the neighbouring magnetic interfaces are oscillating out of 
phase and the Alfven wave becomes " phase- mixed" . This phenomenon decreases dramat- 
ically the perpendicular wave length and generate very sharp transverse gradients. In the 
presence of finite viscosity or resistivity, the short wavelengths are subject to enhanced 
dissipation. Since this mechanism was proposed by Heyvaerts & Priest (1983) for heat- 
ing the solar corona, the linear regime of Alfven wave phase mixing has been intensively 
investigated by many researchers, mainly in the coronal heating context (Abdelatif 1987, 
Browning 1991, Ireland & Priest 1997, see also the critical paper of Parker 1991), includ- 
ing investigation of Alfven wave phase mixing in inhomogencous steady flows (Ryutova & 
Habbal 1995, Nakariakov, Roberts & Murawski 1998) and in two dimensional structures 
(Hood, Ireland & Priest 1997, Ruderman, Nakariakov, & Roberts 1998, Ruderman et al. 
1999, De Moortel et al. 1999, De Moortcl, Hood & Arber 2000) and studies of secondary 
Kelvin-Helmholtz instabilities (Browning & Priest 1984). Possible manifestation of this 
effect in coronal observations have been discussed (Ireland 1996, De Moortel & Hood 
2000) and numerical simulations of Alfven wave phase mixing have been undertaken in 
various astrophysical and geophysical contexts (Botha et al. 2000, Grappin, Leorat & 
Buttighoffer 2000, Malara, Primavera & Vehri 1996, Ofman & Davila 1997). 

In the linear regime of Alfven wave phase mixing, the wave does not interact with 
other wave modes, and the wave dissipation is enhanced due to the generation of small 
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scales across the field. In the nonlinear regime, Alfven waves interact with, in particular, 
magnetoacoustic waves. For the Alfven wave amplitudes expected the lower corona of 
the Sun, the linear treatment of the process of phase mixing seemed to be justified (e.g. 
Sakurai & Granik 1984) as the Alfven waves are subject to cubic, not quadratic as for 
magnetoacoustic waves, nonlinearity. The presence of structuring can dramatically affect 
the efficiency of nonlinear processes. In coronal holes, outwardly propagating Alfven 
waves are amplified by stratification and, therefore, quickly reach the nonlinear regime 
(Ofman & Davila 1997, Nakariakov, Ofman & Arber 2000). In the presence of structuring 
across the field, phase mixing leads to nonlinear generation of transversely propagating 
fast magnetoacoustic waves (Nakariakov, Roberts & Murawski 1997, 1998). 

The nonlinear generation of fast waves, based upon a combination of inhomogeneous 
and nonlinear effects, can provide new physical mechanisms for MHD wave dissipation. 
In the linear regime of phase mixing, Alfven waves dissipate due to the shear component 
of the viscosity tensor, while the fast magnetoacoustic waves dissipate due to the bulk 
viscosity. In the nearly collisionless coronal plasma the bulk viscosity is estimated to be 
10^° times larger than the shear viscosity. Thus, fast magnetoacoustic waves that are 
generated by Alfven waves dissipate much more efficiently than the Alfven waves. This 
mechanism can be regarded as the indirect heating of plasma by phase mixing. 

In the corona of the Sun there is some evidence that the shear viscosity may actually 
be of the same order as the bulk viscosity (Nakariakov et al. 1999). In this case, the 
mechanism of indirect heating is still very important, because it provides a possibily 
of heat distribution across the magnetic field. Indeed, the transversely propagating fast 
waves carry the energy across the field and deposit it at some distance from the region 
of phase mixing. Regular thermal conduction is depressed in the perpendicular direction 
by the magnetic field. Thus, the mechanism of indirect heating of astrophysical plasmas 
by Alfven wave phase mixing requires serious investigation. 

Nakariakov, Roberts & Murawski (1997), by applying a weakly nonlinear description 
for the dynamics of the generated fast magnetoacoustic waves, pointed out that phase 
mixed Alfven waves had to generate secularly growing fast waves (their Eq. (28)). This 
analytical result was successfully reproduced in fuU-MHD numerical simulations of the 
initial stage of the fast wave generation. Botha ct al. (2000) undertook numerical sim- 
ulations of this phenomenon over much longer times, investigating the efficiency of the 
indirect heating by considering harmonic phase- mixed Alfven perturbations. They found 
that initially the amplitude of the generated fast magnetoacoustic waves grows but then 
saturates and does not reach a substantial fraction of the Alfven wave amplitude. They 
further postulated that this saturation was due to the nature of the source term for 
fast waves. For harmonic Alfven waves, as distinct from the single pulse considered in 
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this paper, phase mixing leads to there being many fast wave sources in the direction 
of the background inhomogeneity. Indeed the number of such sources, and their spacing, 
changes in time and as a consequence the fast wave sources cannot remain coherent and 
destructive wave interference hmits the growth of fast wave energy. 

In this work, we simulate numerically the interaction of a linear Alfvcnic pulse with a 
one-dimensional, perpendicular to the magnetic field, plasma inhomogeneity. The pulse 
phasemixes leading to steep transverse gradients. This leads to the generation of fast 
magnetoacoustic waves. This work has three objectives: to investigate whether a phase- 
mixed Alfvcnic pulse, as opposed to Botha et al. 2000's harmonic Alfvcn wave, can 
efficiently generate fast magnetoacoustic waves; to check whether the rapid growth (in 
time) of the generated fast magnetoacoustic waves still holds at the later stages and to 
clarify the role of the destructive wave interference for an isolated pulse. 

The paper is organized as follows: in section 2 we outline the basic aspects of our 
model. Section 3 presents the results of our numerical calculation. While in section 4 we 
close with a brief discussion of the main results of this paper. 

2. The model 

In our model we use equations of cold (zero plasma-/?) ideal MHD 

dV 1 

p—+p{V-\/)V = -—BxcnTlB, (1) 

— = curl(F X B), (2) 

divB = 0, (3) 

^+dW{pV)^0, (4) 

where B is the magnetic field, V is plasma velocity, and p is plasma mass density. 

We consider equations (l)-(4) in Cartesian coordinates {x,y,z) and for simplicity 
assume that there is no variation in the y-direction, i.e. {d/dy = 0). The equilibrium 
state is taken to be an inhomogeneous plasma of density po {x) and a uniform magnetic 
field Bq in the z-direction. 

Nakariakov, Roberts & Murawski 1997 have obtained a set of governing equations 
for the finite amplitude perturbations for the above formulated physical system (see, 
also, Nocera, Priest & HoUweg 1986), and for brevity we do not repeat their derivation 
here. However, we quote their final analytic result, which is two coupled, non-linear wave 
equations describing non-linear coupling of Alfvcn and fast magnetoacoustic waves 



df^ ^' ' dz^ poix) \ dt Air d 
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where Ca{x) ~ Bq/ y^4TrpQ{x) is a local Alfven speed and following notation has been 
used 

9V^ /_ /__N , ^^f.rdV^ , .r9V^\ , B, fdB^ dB,\ BydBy 



^K . .N ^f.rdVy ,,dVy\ B^dBy B^dBy 

A^4 = -^(T4S,-T4B,), (9) 

^5 = -^^(^-^y - ^yBz) + ^("^^y^- - ■^^-^y)' (10) 

Ne^-^{V,B,~V,B,). (11) 

ox 

As noted in Nakariakov, Roberts & Murawski 1997, the non-linear dynamics of Alfven 

waves is described by Eq.(5). Whereas, Eq.(6) is the governing equation for the non-linear 

dynamics of finite amplitude fast magnetoacoustic waves. In the linear limit when the 

right-hand-side terms in the Eqs.(5) and (6) tend to zero, Alfven and fast magnetoacoustic 

waves propagate along the applied magnetic field independently. Therefore, in this case 

Eq.(5) describes phase-mixed Alfven waves, which propagate along 2;-axis perturbing 

physical quantities in the xOz plane. Whereas, Eq.(6) describes either trapped or leaking 

fast magnetoacoustic waves (depending whether there is a maximum or minimum in the 

density profile po{x)). 

We consider the physical situation when initially fast magnetoacoustic perturbations 

are absent and the initial amplitude of the Alfven wave is small. In this case, if the 

Alfven wave is initially a plane wave, the subsequent evolution of the wave, due to the 

difference in local Alfven speed across the x-coordinate, leads to the distortion of the 

wave front. Hence the appearance of transverse (with respect to the applied magnetic 

field) gradients, which grow linearly with time. Therefore the leading non-linear term 

on the right-hand-side of Eq. (6) is the transverse gradient of the magnetic pressure 

perturbation, i.e. to a fairly good accuracy (which is substantiated by our numerical 

calculations presented below) the dynamics of the fast magnetoacoustic waves can be 

described by an approximate equation (Nakariakov, Roberts & Murawski 1997) 

3. Numerical Results 

As in the work of Botha et al. 2000, we have used the following background density profile 

po(x) = 3-2 tanh(Ax) (13). 

Here, A is a free parameter which controls the steepness of the density profile gradient. In 
our normalization, which is the same as that of Botha et al. 2000, the local Alfven speed 
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Fig. 1. Dimcnsionless Alfven speed Ca{x) versus x for A = 0.5 solid line, A = 0.1 long 
dashed line, A = 2.5 short dashed line. 



is Ca{x) — l/-\/3 — 2 tanh(Aa;), and we plot it in Fig.l for the parameters used in our 
numerical simulations (see below). As can be seen from the plot the local Alfven speed 
has a gradient at x = and its value is controlled by the A-parameter. Full numerical 
calculation of the non- linear set of MHD equations (l)-(4), in the above mentioned ge- 
ometry, has been performed using Lare2d (Arber et al. 2001). Lare2d is a numerical code 
which operates by taking a Lagrangian predictor-corrector time step and after each La- 
grangian step all variables are conservatively re-mapped back onto the original Eulerian 
grid using Van Leer gradient limiters. This code was also used to produce the results in 
Botha et al. 2000. 

We set up the code in such a way that initially fast magnetoacoustic perturbations 
are absent and the initial amplitude of the Alfven pulse is small, i.e. A=0.001. In all our 
numerical runs plasma /3 was fixed at 0.001. In the numerical simulations the Alfven per- 
turbation is initially a plane (with respect to x-coordinate) pulse, which has a Gaussian 
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Fig. 2. Solution of the full MHD equations (l)-(4) using Lare2d code, By, at < = 10. 



structure in the z-coordinate, and is moving at the local Alfven speed Ca{x) 

{z -CA{x)tf ' 
5 



By{x, z, t) = Acxp — 



(14). 



Here, (5 is a free parameter which controls the width of the initial Gaussian Alfven pulse. 
The simulation box size is set by the limits —17.5 < a; < 17.5 and —27.5 < z < 27 .h 
and the pulse starts to move from point z — —25.0 towards the positive z's. We have 
performed calculation on various resolutions in order to achieve convergence of the results. 
The graphical results presented here are for the spatial resolution 4400 x 1400, which refers 
to number of grid points in z and x directions respectively. In Fig. 2 we show how the 
initial pulse {By) has evolved hy t — 10. It can be seen from the graph that because of 
the difference in local Alfven speeds (recall from Fig. 1 that for a; < Ca{x) = l/\/5 and 
for x > Ca{x) = 1) the initially flat (with respect to x-coordinate) Gaussian pulse has 
been distorted along x — Q axis. This distortion of the pulse front creates a transverse 
gradient [dBy/dx) which according to Eq. (12) is a driving force for the generation of 
fast magnetoacoustic wave {Vx). At time t — Q fast magnetoacoustic wave was absent 
from the system. While at time t = 10 we see from Fig. 3 that Vx has grown up to the 
square of initial Alfven perturbation amplitude i.e. 10^®. 

Fig. 4 presents a contour plot of By at time t — 50. We observe that the Alfven 
pulse has been distorted even further, thus creating even stronger transverse gradients 
which drive the growth of fast magnetoacoustic waves. In Fig. 5 we plot generated Vx at 
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Fig. 3. Solution of the full MHD equations (l)-(4) using Lare2d code, T4, at i = 10. 

the same time instance, and we observe that its maximal value has grown further up to 
1.78 X 10-6. 

We have kept track of the maximal value of generated fast magnetoacoustic wave, 
^x(t) — max{\Vx{x, z,t)\), as the simulation progresses. The results are presented by 
the think solid line in Fig. 6. As we mentioned above, at the early stages of evolution, 
the dynamics of generated fast magnetoacoustic wave is governed by Eq.(12). In order 
to prove this, we have written a code which, instead of solving the fully non-linear set 
of MHD equations (l)-(4), solves the weakly non-linear Eq.(12) directly, using a second 
order MacCormak scheme, for the driver given by Eq.(14). This driver can be easily 
written by substituting Eq.(14) into Eq.(12): 



z,(.,.,0.-M!ciw''^-<^' 



dx 



exp 



-(z - CAix) tfCAix) t + iz- 2CAix) t) 





2{z-CA{x)tf 



(15) 



where D{x, z, t) is the right hand side of Eq (12). We have solved Eq.(12) on a simulation 
box —30 < a; < 30 and —30 < z < 30 using different spatial resolution. The results for 
the V^{t) are plotted in Fig. 6. We observe that for 4000 x 4000 resolution (thin solid 
line) the behavior of the maximal value of generated fast magnetoacoustic wave is quite 
similar to the results produced by the fully non-linear Lare2d code. In Fig. 7 we plot 
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Fig. 4. Contour plot of solution of the full MHD equations (l)-(4) using Lare2d code, 
By, at t = 50. 

generated Vx at time instance t = 50 obtained by the 2D MacCormak code which solves 
Eq.(12). The similarity with Fig. 5 is obvious, which validates our earlier statement that 
Eq.(12) describes the dynamics of the generated fast magnetoacoustic wave to a sufficient 
accuracy even at later times. 

As can be seen from Eq.(15) the driver of Eq.(12) is a third order polynomial in time. 
Thus, one would expect continuous secular growth of V^it). However, as can be seen from 
Fig. 6, Vx{t) at the later stages grows more slowly than initially. The situation is quite 
similar to what has been observed by Botha et al. 2000. Namely, in spite of secularly 
growing driver the solution of Eq.(12) saturates. Botha ct al. 2000 have speculated that 
the saturation of fast magnetoacoustic wave is due to the destructive wave interference 
effect. Thus, in order to bring clarity into the situation, following Botha et al. 2000, we 
also did calculation of the ID analog of Eq.(12). We have solved equation 

where Di{x,t) = D{x, {Ca{+oo) — Ca{—oo)) t, t). In other words we put z -^ [C a{+oo) — 
Ca{—oo)) t, which physically means that we reduce 2D problem into ID problem by mak- 
ing a transformation to the frame which moves with the mid-point of phase mixed Alfven 
pulse moving with velocity (Cyi(+oo) — C'a{—oo)) t. We solved Eq.(16) with various spa- 
tial resolutions to make sure that the results converge. In Fig. 8 we plot the solution 
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Fig. 5. Solution of the full MHD equations (l)-(4) using Lare2d code, T4, at i = 50. 
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Fig. 6. Maximal value of the generated fast magnetoacoustic wave, V^{t), versus time. 
Thick solid curve represents the solution of the full MHD equations (l)-(4) using Lare2d 
code. Thin solid, long dashed and short dashed curves represent solutions of the Eq.(12) 
using 2D MacCormak code on 4000 x 4000, 1200 x 1200 and 600 x 600 spatial resolution 
accordingly. 
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Fig. 7. Solution of the Eq.(12) using 2D MacCormak code with 4000 x 4000 spatial 
resolution, T4, at t = 50. 



of Eq.(16) on the simulation interval —55 < a; < 55 using 20000 grid points. At time 
t = t4 is zero. Top row shows Vx at time instances t — 1,10,50. We see that the 
solution comprises two pulses outgoing into positive and negative directions at speeds 
Ca{x) = l/\/5 for X < and Ca{x) = 1 for x > 0. The bottom row contains plots of 
the driver, Di{x,t), at the same instances. Interestingly, we observe that initially, the 
driver is asymmetric and has larger positive part, however at the later stages it becomes 
symmetric and progressively narrow. This means that, similarly to results obtained by 
Botha et al. 2000, the fast magnetoacoustic waves generated by such a driver will saturate 
because the positive and negative parts of the driver, which get closer as time progresses, 
cancell each other due to destructive wave interference. The difference between our re- 
sults and the ones by Botha at al. 2000 is that since we have considered a single Gaussian 
pulse (hoping that single pulse will not trigger destructive wave interference effect) we 
have a single symmetric driver. Botha at al. 2000 considered a harmonic initial Alfvcn 
perturbation, thus obtained a driver with a number of symmetric peaks, who's number 
grows in time (see Fig. 6 in Botha et al. 2000 for comparison). We have kept track of the 
maximal value of generated fast magnetoacoustic wave, V^it), as our ID code progressed 
in time. The results for the three values of A are presented in Fig. 9. It should be noted 
that there is a notable difference between 2D (Fig. 6) and ID (Fig. 9) results. In 2D, 
as opposed to ID, we do not observe clear saturation of V^it). We can only conclude 
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Fig. 8. Top row: solution of Eq.(16) using ID MacCormak code on 20000 spatial res- 
olution, Vx, at times t = 1,10,50; bottom row: driver of Eq.(16), Di{x,t), at times 
t= 1,10,50. 



form Fig. 6 that at the later stages V^{t) grows much slower than initially. This can 
be explained by the fact that in ID, we define V^{t) as a maximal value of generated 
fast magnetoacoustic wave in the frame which moves with the mid-point of phase mixed 
Alfvcn pulse having the velocity (C^(-l-oo) — C^(— oo))t. Therefore in 2D the quantity 
that would be more appropriate to compare to V^{t) from ID results is the maximal 
value of generated fast magnetoacoustic wave in ID slices along x-coordinate which are 
taken at z — (Ca(+oo) — Ca(— oo))t as the 2D MacCormak code progress in time. We 
denote the latter quantity by V^{t) =- maa;(|T4(a;, z = (Ca(+oo) — C^(— oo)) i,i)|). In 
Fig. 10 we plot y^(t) for three values of A = 0.1,0.5,2.5. These are the same values 
that are used in Fig.l. The A-parameter controls strength of the transversal gradient. 
The noteworthy features that can be observed from the plot are as following: stronger 
transverse gradients (larger A's) cause earlier start of the fast magnetoacoustic waves 
saturation process, and in turn, yield lower saturation levels. Whereas, weaker gradi- 
ents (smaller A's) result in later onset of the saturation, and in turn, higher saturation 
levels. This result, actually, validates our explanation of the fast magnetoacoustic wave 
saturation by the destructive wave interference effect. Namely, larger A means steeper 
transverse gradient, i.e. narrower driver acting from the beginning. Therefore, destruc- 
tive wave interference effect starts to a play role earlier yielding faster shut down in the 
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Fig. 9. Maximal value of the generated fast magnetoacoustic wave, V^{t), versus time. 
Solid, long dashed and short dashed curves represent solutions of Eq.(16) using ID Mac- 
Cormak code on 200000 spatial resolution for A = 0.5, 0.1, 2.5 accordingly. 



growth of Vx- For the the weaker gradients, i.e. small A, the opposite is true. This result 
is in agreements with that obtained by Botha et al. 2000 (see Fig. 3 in Botha et al. 2000 
for comparison). Another, important observation from Fig. 10 is that the saturation of 
V^{t) is similar to that from ID results (Fig. 9). This differs from the results shown in 
Fig. 6 in which the amplitude of Vx continues to grow albeit not at the rapid initial phase 
from t = to t = 5. The discrepancy between the ID and 2D results can be understood 
from the form of the 2D solution in Fig. 5. Here it is clear that the maximum Vx in 
the box is not on the x line used to produce Fig. 10 but is located at regions where the 
background density levels off to a constant value. At these locations the x gradient of By 
remains small but does not contract in x as rapidly as in the centre of the phasemixing 
region. Thus these regions are a source of continued secular growth in Vx for longer. 

Following Botha et al. 2000 we have performed a study of parameter space, by inves- 
tigating the dependence of the fast magnetoacoustic wave saturation levels upon the two 
free parameters in our problem, namely, A and 6, using the 2D MacCormak code. In Fig. 
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Fig. 10. Maximal value of the generated fast magnetoacoustic wave, V^{t), in ID slices 
along x-coordinate which are taken at z — {Ca{+oo) — C/i(— oo))t versus time. Solid, 
long dashed and short dashed curves represent solutions of Eq.(12) using 2D MacCormak 
code on 4000 x 4000 spatial resolution for A = 0.5, 0.1, 2.5 accordingly. 



11 we plot V^{t) for three values of A = 0.1, 0.5, 2.5. Again, the noteworthy features are: 
stronger transverse gradients (larger A's) cause earlier start of the fast magnetoacoustic 
waves saturation process, and in turn, yield lower saturation levels. Whereas, weaker gra- 
dients (smaller A's) result in later onset of the saturation, and in turn, higher saturation 
levels. It is worthwhile to note the ID results (Fig. 9) of the variation of A-parameter 
show a similar trend, i.e. start of the saturation process earlier for larger A's. However, 
as in the ID results of Botha et al. 2000 (sec p. 1191 in their paper for the discussion) 
the scaling of the saturation amplitudes with A is not the same as in 2D or full MHD 
simulations. 

In Fig. 12 we investigate the scaling of saturation levels of the generated fast magne- 
toacoustic waves with the (5-parameter, which controls width of the phase-mixed Alfven 
pulse over z-coordinate. Fig. 12 confirms that with the decrease of S, i.e. for narrower 
Gaussian pulses, saturation levels of the generated fast magnetoacoustic waves increase 
and vice versa. 

Finally, we would like to comment on the choice of the upper bound of our simulation 
time, which we have set to t =: 50.0. DcMoortel et al. 1999 have found that in linear 
phase mixing (with the dissipation effects included) shear viscosity dissipation reaches 
its maximum when the harmonic Alfven wave is 10 wavelengths out of phase. Assuming 
that this result holds for the traveling Gaussian Alfven pulse, we can estimate time the 
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Fig. 11. Maximal value of the generated fast magnetoacoustic wave, V!^(t), versus time. 
Solid, long dashed and short dashed curves represent solutions of Eq.(12) using 2D Mac- 
Cormak code on 4000 x 4000 spatial resolution for A = 0.5, 0.1, 2.5 accordingly. 
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Fig. 12. Maximal value of the generated fast magnetoacoustic wave, V^{i), versus time. 
Solid, long dashed and short dashed curves represent solutions of Eq.(12) using 2D Mac- 
Cormak code on 4000 x 4000 spatial resolution for 6 = 0.1,0.005,0.5 accordingly. 



i» it takes for the classic shear viscosity to become significant. The distance between the 
wings of the distorted phase mixed Alfven pulse (see Fig. 4) is (C^(+oo) — C/i(— oo)) t*. 
If we use the pulse width, 2\/(5, instead of wavelength in the original DeMoortel et al. 
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1999's estimate, we can write (Cyi(+oo) — C/i(— oo)) t* = 10 x (2\M). Inserting our values 
lor the asymptotic Allvcn speeds and 5 = 0.1 in the latter equality, we deduce t* « 11. 
Therefore, there is little point running the present model for t > i*, because as shown 
by DeMoortel et al. 1999, viscosity would be significant for such later times and this is 
excluded from our current model. 

4. Discussion 

The motivation ol this paper was three-fold. Botha et al. 2000 concluded that for har- 
monic phase-mixed Alfvcn perturbations the generated fast magnetoacoustic waves do 
not grow to a substantial fraction of the Alfvcn wave amplitude. Firstly, we wanted to 
test whether a single phase-mixed Gaussian Alfvcn pulse behaves in the same way. Sec- 
ondly, to test the validity of the weakly non-linear analytic description of Nakariakov, 
Roberts & Murawski 1997 for large times. Thirdly, to check whether the rapid growth 
(in time) of the generated fast magnetoacoustic waves still holds at the later stages. 
Our present analysis has clearly demonstrated that: 

— Small-amplitude, phase-mixed, single Alfvcn pulses do not generate fast magnetoa- 
coustic waves with large amplitudes. The saturation is due to the destructive wave 
interference. 

— We have confirmed that the weakly non-linear analytic description of Nakariakov, 
Roberts & Murawski 1997 (Eq.(12)) is valid and produces results that are entirely 
consistent with the simulation results from fully non- linear MHD equations (l)-(4) at 
all stages, provided the initial amplitude of the Alfvcn perturbation is small. 

— We have found that the fast magnetoacoustic wave, generated by the phase mixed 
single Alfvcn pulse, at the later stages grows much slower than at earlier times due 
to the destructive wave interference. 

We have also established that the stronger transverse gradients (larger A's) cause 
earlier start of the fast magnetoacoustic waves saturation process, and in turn, yield 
lower saturation levels. Whereas, weaker gradients (smaller A's) result in later onset 
of the saturation, and in turn, higher saturation levels. This result, actually, validates 
our explanation of the fast magnetoacoustic wave saturation by the destructive wave 
interference effect. Namely, larger A means steeper transverse gradient, i.e. narrower 
driver, acting from the beginning. Therefore, destructive wave interference effect starts 
to a play role earlier yielding faster shut down in the growth of Vx. For the the weaker 
gradients, i.e. small A exactly the opposite is true. This result is in agreements with 
the one obtained by Botha et al. 2000 for a harmonic Alfvcn wave (see Fig. 3 in Botha 
et al. 2000 for comparison). It is worthwhile to mention, however, that on one hand, 
contrary to the one's expectations, steeper inhomogenities of the Alfven speed do not 
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provide the higher saturation levels of the fast wave generation. In other words, with 
increase of A saturation levels decrease. On the other hand, in the homogeneous plasma 
(A = 0) there is no generation of the fast magnetoacoustic wave by a plane (with respect 
to x-coordinate) Alfven wave at all, i.e. V^{t) -^ when A ^ 0. Therefore, there is a 
certain optimal value of gradient of the inhomogcneity for which the fast wave saturation 
level is maximal. At this stage we have not performed full parametric study in order to 
determine the optimal value of A due to the vast requirements for CPU time. 

The mechanism of indirect plasma heating by the fast magnetoacoustic waves that are 
non-linearly generated by the phase-mixed Alfven perturbation, discussed in this paper, 
is certainly applicable for the coronal heating problem through dissipation of fast mag- 
netoacoustic waves in open (e.g. plumes) and closed (loops) coronal plasma structures. 
Although, we found that the proposed mechanism, in the considered set up, is limited 
by destructive wave interference, further work is needed to find possible ways around of 
the encountered difficulty. Possible cures might be considering 2D density inhomogcneity 
profiles, that would prevent the driver of Eq.(12) from becoming too narrow, and/or in- 
clusion of the dissipation into our model, that would prevent the transverse scale collapse, 
which all ultimately will lead to the decrease of the destructive wave interference. 

Acknowledgements. Numerical calculations of this work were done using the PPARC funded 
Compaq MHD Cluster in St. Andrews. 
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